Do multineutrons exist? 
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Following recent work in which events which may correspond to a bound tetraneutron ( 4 n) were 
observed, it is pointed out that from the theoretical perspective the two-body nucleon-nucleon (NN) 
force cannot by itself bind four neutrons, even if it could bind a dineutron. Unrealistic modifications 
of the NN force or introduction of unreaslistic four-nucleon force would be needed in order to bind 
the tetraneutron. The existence of other multineutron systems is discussed. 

PACS numbers: 21.10.Dr, 27.10.+h, 27.20.+n, 
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In a recently reported experiment [l| events were observed that exhibit the characteristics of a multineutron cluster 
C*~) 1 liberated in the breakup of 14 Be, most probably in the channel 10 Be+ 4 n. The lifetime of order 100 ns or longer 
suggested by this measurement, would indicate that the tetraneutron is particle stable. The existence of a bound 
tetraneutron, if confirmed, could challenge our understanding of nuclear few-body systems and nucleon-nucleon (NN) 
interactions. In the present letter, I would like to make several observations in this context. 

a: i 

The tetraneutron-like events seen in [| deserve attention because the breakup 14 Be — > 10 Be+ 4 n represents one of 
the best possible tools to search for a tetraneutron. In earlier experiments, the tetraneutron was searched for using 
heavy-ion transfer reactions such as 7 Li( n B, 14 0) 4 n 2], 7 Li( 7 Li, 10 C) 4 n 0] and double exchange reaction 4 He(7r~ , 7r + ) 4 n 
Q, |3- These reactions require considerable reconfiguring of the target nuclei and should be strongly suppressed. 
A negative outcome from these reactions could be easily anticipated. On contrary, the nucleus 14 Be consists of a 
strongly bound core 10 Be and four valence neutrons whose separation energy is only about 5 MeV. In the attractive 
field of 10 Be these four neutrons could form a tetraneutron-like configuration which might be shaken off in the 14 Be 
breakup. 
II 

No proper ab-initio four-body calculations of the tetraneutron with realistic two-body and three-body NN forces 
are known to the author. However, several calculations of the tetraneutron are available. 

(i) In Ref. [| the tetraneutron, studied in the translationally invariant AHlo oscillator shell model, was found to be 
unbound by 18.5 MeV. However, the oscillator basis is not appropriate for the description of unbound systems. 

(ii) No bound tetraneutron was found in Ref. within the angular potential functions method with semirealistic 
Q 1 NN interactions. No search for the four-body resonance state was carried out. 

(Hi) No bound tetraneutron was found in Ref. within the stochastic variational method on a correlated Gaussian 
basis for a range of simple effective NN potentials. No search for the four-body resonance state was carried out. 
J> \ (iv) A search for four-body resonances in the lowest order of the hyperspherical functions method (HSFM) gave a 
•'-j ■ null result 0, 

(v) The energy behaviour of the eigenphases, studied in Ref. within the HSFM using the K max = 6 model 
space, led the authors to the conclusion that the tetraneutron may exist as a resonance in the four-body continuum 
at an energy between 1 and 3 MeV. However, a clear indication of the resonance has been seen with only one of 
the NN potentials used in the calculations, namely with the Volkov effective NN force VI Volkov effective NN 
interactions reproduce the experimental binding energy of another four-body system, 4 He, but they bind a dineutron. 
In the particular case of VI, the dineutron is bound by 0.547 MeV. Therefore, VI cannot be used in calculations 
of multineutron systems. Another NN potential, used in Ref. 11], namely that of Reichstein and Tang (RT) [l3j . 
reproduces the n-p triplet and p-p singlet scattering lengths and does not bind a dineutron. With this potential the 
energy derivatives of the eigenphases monotonically decrease with increasing energy. This means that resonances in the 
four neutron system are absent, at least within the model space considered. Since no convergence of the eigenphases 
with an increase of the hyperangular momentum has been achieved in these calculations, such a conclusion does not 
look convincing. 

To understand whether the RT potential can produce any resonance or bound state if the model space of the HSFM 
is increased, I have studied the tetraneutron within the extreme uncoupled adiabatic approximation of the HSFM 
[l4| which provides a lower limit for the binding energy ^| within the model space considered. The model space in 
these calculations has been increased up to K max = 16. In this limit it is sufficient to diagonalise the matrix of the 
hyperradial potentials Vx'-y 1 ,Ky{p), which include both nuclear and centrifugal force, and to solve a Schrodinger-like 
equation with the lowest diagonalised potential Vdi ag (p) E3- The hyperradial potentials have been calculated using 
the technique developed in Ref. [l6| . 

The RT potential has only central triplet and singlet even components Vte and Vse and the singlet and triplet 
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FIG. 1: (a) The lowest diagonalised hyperradial potentials Vdiag(p) of 4 n, 6 n, 8 n and 10 n calculated with Kmax = 2, 4, 16, 
K max = 4, 6,..., 12, A" m ai = 6,..., 12 and K max = 10, 12, 14 respectively for the RT potential; (6) the same as in (a) only for 
4 n and 6 n and the VI potential. 



odd potentials are usually obtained from them using an arbitrary parameter u so that Vso = {u — 1)Vse and 
Vto = (u — 1)Vte- The present calculations have been performed with u = 1 which corresponds to no interaction 
in the odd partial waves. These calculations provide a lower limit for the binding energy because realistic values of 
u, which lie in the range of 0.93 to 0.98, induce repulsion between some of neutrons, which can only decrease the 
binding. 

The diagonalised hyperradial potentials Vdiag(p), calculated for different model spaces K max , are plotted in Fig. la 
as a function of hyperradius p. They have almost converged and are purely repulsive, monotonically decreasing 
with the hyperradius, without any sign of local attractive pockets. Therefore, the RT potential can neither bind the 
tctrancutron nor produce any resonances. To get a bound tetraneutron, the value of u should be increased to u = 2.3. 
Such a value corresponds to a unphysical attraction in odd partial waves and binds the isotopes 4 H and 5 He with 
respect to the t + n and 4 Hc + n thresholds by about 6.5 MeV and 32 MeV respectively. In reality, these nuclei are 
unbound by 3 and 0.9 MeV. 

For comparison, the calculations of the hyperradial potentials have been performed with Volkov NN potential VI. 
The diagonalised hyperradial potentials Vdi ag (p), shown in Fig. lb, reveal local attractive pockets for K max > 6. 
These pockets become negative for K max > 12, but they are too shallow to form a bound state. To get a bound 
tctrancutron, the Majorana parameter m should be changed from its standard value of 0.6 to m = —0.2. Such a 
change provides E( 4 n) = —1.2 MeV and does not influence the binding energy of 4 He, however, it binds 4 H and 5 He 
with respect to the t + n and 4 He + n thresholds by 7 MeV and 25 MeV respectively. 

Ill 

The results of the theoretical calculations suggest that the tetraneutron must be unbound and most likely should 
not exist as a resonance, at least when only two-body central forces are considered. It is remarkable that even when an 
effective NN potential binds the dineutron (the case of VI), it still cannot bind two dineutrons, although a resonance 
in such a system should exist. In order for the two-body central force to be able to bind the tetraneutron, a huge 
unphysical attraction must be introduced in the triplet odd potential. Such an attraction would strongly overbind 
4 H, 5 He and other known A > 4 nuclei. The necessity for such a huge attraction lies in the relative number of nucleon 
pairs in even and odd states as determined by the Pauli principle. The NN force binds two nuclcons only in the triplet 
even state. It is not sufficiently attractive in singlet even states, and is repulsive in the triplet and singlet odd states. 
The probabilities P to find a pair of nucleons in all these states are shown in Table 1 for the dominant components 
[22] 51 S', [31] 33 P and [4] n 5 of 4 n, 4 H and 4 He respectively, together with binding energies of 4 H and 4 He. One can 
see that the binding energy of a nucleus is strongly correlated with the probabilities to find a pair of nucleons in even 
states, especially in the triplet even state where the n — p bound state exists. In 4 H, a decrease of the probabilities 
Pte and Pse by 30% in comparison with 4 He leads to a dramatic decrease of its binding energy by about 20 MeV. 
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TABLE I: Binding energies (in MeV) and the probabilities P to find a pair of nucleon in the triplet even (TE), singlet even 
(SE), singlet odd (SO) and triplet odd (TO) states for different A — 4 isotopes and for 6 ' 8 ' 10 n. 
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In the tetraneutron, triplet even components (the only ones which produce binding in the NN system) are absent 
and the probability to find two neutrons in the repulsive triplet odd state increases by 50% in comparison with 4 H. 
Therefore, it is not unreasonable to assume that moving from 4 H to 4 n one can loose a comparable amount of binding 
energy as going from 4 Hc to 4 H. This would definitely make the tetraneutron unbound. 
IV 

The conclusion of particle instability of 4 n made above has been drawn on the basis of only two-body NN forces. 
The presence of the 3N and, possibly, 4N forces could, in principle, modify this conclusion. The contribution of the 
3N potential to the binding energies of neutron drops n and 8 n formed in an external field has been calculated with 
Green's function Monte Carlo Methods in Ref. 0]. It ranges from 1 to 5 MeV. A similar contribution from the 3N 
force could be expected in the case of 4 n. But such a contribution will not significantly lower the hyperradial potential 
from Fig. la to produce a potential well in which a bound state of four neutrons could be formed. 

As for the 4N force, the results of the ab-initio calculations of the 3 < A < 8 nuclei suggest that either the 4N 
contributions to the observed binding energies of these nuclei are smaller than 1%, or their effects are included in some 
parts of the 3N forces 0J ■ The latest models of the 3N force used to calculated the binding energies of these nuclei 
have not yet been tested in the description of polarization observables of the low-energy Nd scattering where the 
contribution of the 3N force is important 0], [13. Besides, the latest numerically accurate calculations employing a 
3N force do not reproduce the proton analyzing power for the p- 3 He scattering |20j . The simultaneous fit of the Nd 
and p- 3 He polarization observables and of binding energies of the lightest nuclei may leave room for a 4N force. 

In order to get an idea of how strong the 4N force should be to bind the tetraneutron, the HSFM calculations of the 
present paper have been repeated with the RT potential using a 4N force simulated by the potential V~4n(p) — Woe~ ap . 
The values a = 0.7, 1.2 and 1.5 /to -1 used in these calculations, were the same as the range of the phenomenological 
term of the 3N spin-orbit force introduced in Ref. . The tetraneutron becomes bound if the corresponding values 
Wo are equal to —410, —1460 and —2530 MeV respectively. These values are two orders of magnitude larger than the 
values —1, —10, and —20 MeV obtained for the 3N spin-orbit force (lfij. For the Volkov potential VI, these strengths 
are similar: Wq = —320, —1565 and —2900 MeV. If the same 4N force existed in 4 He, the binding energy of 4 He, 
calculated with VI, would be —88, —82 and —134 MeV respectively instead of —28.3 MeV, which would be impossible 
not to notice. 

V 

The reported observation of the tetraneutron revives an old question "Do multineutrons exist?" Theoretical inves- 
tigations of multineutrons have been carried out only for 6 n and 8 n. The 6hu> oscillator shell model, with the same 
effective interactions that unbind the tetraneutron by 18.5 MeV, predicts that 6 n is unbound only by 6.5 MeV [2l| . 
However, the with the same interaction, the hydrogen isotope 6 H should be unbound with respect to the t + n + n + 
n decay only by 1.3 MeV and thus should be less unbound than recently discovered isotope 5 H j2^. It this was true, 
the 6 H isotope would have been already observed. Other calculations of multineutrons, performed within the angular 
potential functions method 0, predict that 6 n together with another multineutron 8 n are unbound. No search for 
six- or eight-body resonances have been carried out in these works. 

In the present letter, the hyperradial potentials Vdiag{p) have been calculated for 6 n with the same NN potentials as 
for 4 n, for K max ranging from 4 to 12. The calculations with the RT potential and u = 1, shown at Fig. la, imply that 
6 n should be even more unbound than n and that it can not exist as a six-body resonance. The reason for this lies 
in further decrease of the number of the singlet even NN pairs, where attraction occurs, and increase of the number 
of repulsive triplet odd pairs (see Table 1). In the case of the NN potential VI, which binds the dineutron, the model 
space is not large enough to make definite conclusions about existence of 6 n. It is possible that further increase of the 
model space will create a local attraction pocket where a six-body resonance could be formed but it is unlikely that 
it will lead to a six neutron bound state. If repulsion in the triplet odd components was replaced by attraction in 
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such a way that 4 n is bound (u — 2.3 case for the RT potential and m — —0.2 case for VI) then 6 n would be bound 
by about 22 MeV or 50 MeV respectively. This results does not look surprising because there is more repulsion in 6 n 
than in n and substitution of the repulsion by attraction should lead to more binding in 6 n than in 4 n. 

Finally, the hypcrradial potentials for 8 n and 10 n have been calculated with the RT potential up to K max equal to 
12 and 14 respectively. These potentials, shown in Fig. la, do not exhibit any sign of either a bound or a resonant 
state which is consistent with further decrease of attraction and increase of repulsion in these systems demonstrated 
in Table 1. 

Summarizing, due to the small probability for a pair of neutrons to be in the singlet even state, the two-body 
NN force cannot by itself bind four neutrons, even if it could bind a dineutron. Unrealistic modifications of the NN 
force would be needed to bind the tetraneutron. Therefore, it is unlikely that the events in the breakup of 14 Be were 
caused by a formation of a bound tetraneutron. A different explanation should be sought for this experiment and 
new experiments arc needed to clarify this issue. As for experimental searches of other multineutron systems, the 
calculations presented here suggest that they might be unsuccessful. 
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